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1 Introduction 

The electroweak penguin mediated B — ► X s £ + £~ decay is a good probe to study be- 
yond the Standard Model physics at low energy The effect of new physics could be 
observed as corrections to the branching fraction, di-lepton invariant mass spectrum, 
and to the Wilson coefficients such as Cg, Cxo, and a phase of C-j [1-10]. Experimen- 
tally, the exclusive B — > K£ + £~ decay was observed by the Belle collaboration for the 
first time [IT] and has also recently confirmed by the BABAR collaboration |T2| . We 



present here a preliminary result on a search for the inclusive B — > X s £ + £ decay, 
which can be calculated more reliably than exclusive decay modes. 



2 Analysis 

The data sample used in this analysis consists of 43 fb _1 of e + e _ collisions at the 
T(AS) resonance collected by the Belle detector [[H| at the KEKB storage ring. 

The signal is identified as X s £ + £~, where X s is defined to be a K ± or K$ accom- 
panying to 4 pions including one n° at most, and t + t~ is an oppositely charged 
electron or muon pair. Charged tracks are required to originate from the interaction 
point and be well identified by the Belle particle identification devices. Electrons 
are identified using the ratio between the energy deposited in the electromagnetic 
calorimeter and the measured track momentum(i?/p). Electrons are required to have 
momenta above 0.5 GeV/c. Muon candidates are required to reach the iron flux re- 
turn and have momenta greater than 1.0 GeV/c. The muon identification efficiency 
is about 80 to 85 % with a fake rate of 1 to 2 %. Charged hadrons are identified by 
the combined response of aerogel counters, time of flight, and dE/dx measured in the 
drift chamber. For charged kaons, we have an efficiency of 85 to 90% and a pion fake 
rate that is less than 10%. Candidate 7r°'s are identified from pairs of photons with 
energy deposits of at least 50 MeV and an invariant mass within 10 MeV/c 2 of the 
nominal 7r° mass. 



Events with leptons from charmonia, such as J/ ip and ip' are vetoed if the dilepton 
mass lies around the nominal charmonia masses. 

Background from qq continuum events is suppressed by requiring the ratio of the 
second to zeroth Fox- Wolfram moments R 2 = H 2 /H < 0.35 and | cos9 t hrust\ < 0.85, 
where Othrust is the angle between the thrust axis of the B candidate and that of the 
rest of the event. The | cos Othrust | distribution of continuum events peaks at 1 while 
that of signal events is flat. 

To reduce backgrounds further and select the best candidate in an event, we 
introduce four kinematic variables: the angles between the K and £ + (0ki+), and the K 
and £~{0 Ke -), the B flight direction^), and the energy difference (AE = E B —E heam ). 
If the K meson and £ + £~ pair originated from a photon or Z boson; since these are 
emitted in opposite directions, the sum of cos 9ki+ and cosOxe- should be negative 
for the signal. The polar angle of B flight direction has a 1 — cos 2 9b distribution 
while background is flat. The signal AE distribution peaks at zero. 

A likelihood ratio LR is determined by parameterizing distributions of three kine- 
matic variables with Monte Carlo simulation for both the signal and background. It 
is expressed as LR = p S ig/(p S ig + Pbg), where p S i g and pbg are probability density 
functions for the signal and background distributions, respectively. An LR selection 
requirement is chosen to maximize the figure of merit S/ V S + N in the signal region 
where S is the expected signal yield assuming the SM prediction and N is the number 
of expected background events. 

The invariant mass of the X s must satisfy Mx s < 2.1 GeV/c 2 in order to reject 
combinatorial background. 

After the application of all selection requirements, the reconstruction efficiencies 
are estimated by Monte Carlo simulation based on the inclusive model in the recent 
paper by Ali et ai.ffl and based on the exclusive B — > K£ + £~ and B — > K*£ + £~ 
model by Greub et al. |2j. The nominal reconstruction efficiencies are determined to 
be 3.6 % for B X s e+e- and 3.8 % for B -> 

Most of the background is combinatorial background from B and D decays. The 
background is determined by fitting the distribution of beam constrained mass M& c . 
The beam constrained mass is expressed as M bc = (Eb ea m) 2 ~ (pb) 2 , where ps is 
the B candidate's center-of-mass momentum vector and Eb eam is the center-of-mass 
energy of B meson, respectively. In this distribution, the signal peaks at the B 
meson mass and is fitted by a Gaussian function while the combinatorial background 
is parameterized by a phase space function with a kinematic threshold (the ARGUS 
function). 

Background from B — > X s 7t + it~ decay can make a small contribution to the signal 
peak in the M bc distribution when both the ir + and ir~ are misidentified as muons. The 
muon misidentification rate is determined to be 1 to 2 % in the laboratory momentum 
range p^ > l.OGeV/c. To estimate the contamination in the signal region, a B — » 
X s n + Tr~ sample is selected by applying all the signal selection requirements except 
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for the muon identification. The yield is multiplied by the momentum dependent fake 
rate for each pion. This B — > X s tt + 7t~ background is estimated to be 2.4 events and 
subtracted from the yield of the signal fit. 

3 Results and discussion 

Figure [L] shows the Mb c distribution for B — > X s e + e~ , B — > X s fi + fi~, and for the 
combined B — > X s £ + £~ samples, respectively. The B — > X s e ± /i T sample is used 
to check the background parameterization. The fit results, signal yields, branching 
fractions, and statistical significances are shown in Table [I]. 

Systematic uncertainties in this measurement are listed in Table [| The largest 
contribution to the systematic uncertainty is the decay modeling for B — ► X s e + e~ 
and B — > X s fi + fi~ decays. 

We measure the inclusive B — > X s £ + £~ decay branching fraction for the first time 
with a statistical significance of greater than 4cr. Preliminary results for branch- 
ing fractions are B(B -> X s ^+fi~) = (8.9j£?l};?) x 10" 6 and B(B -> X s £ + £-) = 
(7.liJ|gii; 2 ) x where the first error is statistical and the second is systematic. 

The branching fractions and upper limits determined by this study agree well with 
theoretical predictions in the SM framework. We expect that both the theoretical and 
experimental errors will decrease so that theories based on SM can be tested with 
precision measurements in the near future. 

The Mx a and dilepton invariant mass plots are shown in Figures |2| and |3|, re- 
spectively. The Mx s distribution extends beyond the K and K* mass region. The 
dips around 3.1GeV/c 2 and 3.7GeV/c 2 in the Mu distribution are due to J/ip and 
ip' vetoes. The experimental and theoretical uncertainties are expected to improve as 
more data is accumulated. Precision measurements can be done for SM and beyond 
SM predictions. 

In summary, evidence for the inclusive B —>■ X s £ + £~ decay has been presented 
for the first time. The experimental and theoretical uncertainties are expected to 
decrease with time. Precision measurements can be done in near future that test the 
Standard Model and Beyond the Standard Model predictions. 
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Mode 


#signal 


B.F.(xlO- 6 ) 


Signif. 


X s e + e~ 
X s £ + £- 


1D - D -7.3-3.8 
on 7+7.9+5.4 
OU. < _7.4_3.8 
a 7 R+ll.Q+9.6 
^' -°-10.4-8.0 


< 11.0 

o q+2. 3+1.6 
°- y -2.1-1.7 
7 -,+1.6+1.4 
' - 1 - 1.6-1.2 


2.1 
4.4 
4.8 



Table 1: Preliminary results for I? — > X s e + e~, 
X s fi + fi~, and combined Significance is 

extracted from the statistical error only. Only a 
90% confidence level upper limit is shown for the 
B — > X s e + e~ process. 



oource 


A s e e 


v ,,+ ,,— 

X s flL^H 


Tracking 


8.1 % 


8.0 % 


ivaoii iu 


1 Q % 


z.u /o 


Pion ID 


8 % 


8 % 


Lepton ID 


3.6 % 


4.4 % 


Ks detection 


2.1 % 


1.5 % 


7T detection 


2.0 % 


1.6 % 


MC stat. 


3.9 % 


4.1 % 


Decay model 


+14% 
-9% 


+16% 
-12% 


Total 


+18% 
-14% 


+19% 
-16% 



Table 2: Systematic uncertainty 
summary. 
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Figure 1: M bc distributions and fit results. Top-left: X s e + e~ candidates, top- right: 
X s /i + fi~ candidates, bottom- left: X s £ + £~ = (X s e + e~) + (X s fi + fi~) candidates, and 
bottom-right: X s e ± fi T to estimate combinatorial background. The significance is 
determined from the statistical error only. 
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Figure 2: Comparison of Figure 3: Comparison of 

the Mx 3 distribution for the Mu distribution for 

data and MC. data and MC. 
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